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Abstract 

Nonlocality of two-mode states of light is addressed by means of CHSH inequality based on 
displaced on/off photodetection. Effects due to non-unit quantum efficiency and nonzero dark 
counts are taken into account. Nonlocality of both balanced and unbalanced superpositions of 
few photon-number states, as well as that of multiphoton twin beams, is investigated. We find 
that unbalanced superpositions show larger nonlocality than balanced one when noise affects the 
photodetection process. De-Gaussification by means of (inconclusive) photon subtraction is shown 
to enhance nonlocality of twin beams in the low energy regime. We also show that when the 
measurement is described by a POVM, rather than a set of projectors, the maximum achievable 
value of the Bell parameter in the CHSH inequality is decreased, and is no longer given by the 
Cirel'son bound. 
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I. INTRODUCTION 



Quantum entanglement for both discrete and continuous variable systems has been exten- 
sively analyzed, also revealing its subtle relations with other quantum mechanical features 
such as nonlocality. Indeed, it has been pointed out that the concept of entanglement coin- 
cides with nonlocality only for the simple case of bipartite pure states. As soon as we deal 
with mixed states, entangled states can be found which do not show properties of nonlo- 
cality while, not unexpectedly, the converse is always true In addition, the amount of 
nonlocality, i.e., the amount of violation of a suitable Bell inequality, crucially depends on 
the nonlocality test adopted in the analysis, ranging from no violation to maximal violation 
for the same (entangled) quantum state. 

In this paper we address nonlocality of different kinds of two-mode states of light by 
means of displaced on/off photodetection taking into account the effects of non-unit quan- 
tum efficiency and dark counts. This kind of measurement was first proposed in Ref. 
where, in particular, it was pointed out that the correlation functions violating the Bell 
inequalities (in the ideal case) involve the joint two-mode Q-function. The reason to pay 
a particular attention to on/off tests of nonlocality is twofold. On one hand, it has been 
shown that violation of Bell inequalities may be quite pronounced for some relevant state 
of light in the ideal case On the other hand, and more importantly, on/off tests may be 
effectively implemented with currently technology. In this framework, it is of interest to take 
into account the effects of experimental imperfections, e.g., non-unit quantum efficiency and 
nonzero dark counts j^, and to investigate the nonlocality properties of physically realiz- 
able entangled states. Indeed, realistic implementations of quantum information protocols 
require the investigation of nonlocality properties of quantum states in a noisy environment. 
In particular, the robustness of nonlocality should be addressed, as well as the design of 
protocols to preserve and possibly enhance nonlocality in the presence of noise. 

The paper is structured as follows: in the next Section we describe in some detail the 
nonlocality test we use throughout the paper, while in Section ITTTl we analyze nonlocality of 
superpositions, both balanced (Bell states) and unbalanced, involving zero- and one-photon 
states. In Section HXl we address the nonlocality of superpositions containing two-photon 
states, whereas Section El is focused on multiphoton twin-beam state. In Section ED we 
analyze the effect of dark counts on the violation of CHSH inequality, whereas in Section IVTTl 
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we address inconclusive photon subtraction (IPS) as a method to enhance nonlocality of twin- 
beam. Section IVIIII is devoted to a more detailed analysis of the choice of parametrization 
leading to violation of inequalities, whereas, in Section |lXl we show how a nonlocality test 
based on POVM measurement cannot yield the maximal violation of inequalities expressed 
by the Cirel'son bound. Finally, Section |X| closes the paper with some concluding remarks. 



II. CORRELATION FUNCTIONS AND BELL PARAMETER 

The nonlocality test we are going to analyze is schematically depicted in Fig.HJ two modes 
of the radiation field, a and b, are excited in a given (entangled) two-mode state described 
by the density matrix g, and then are locally displaced by an amount a and f3 respectively. 
Finally, the two modes are revealed by on/off photodetectors, i.e., detectors which have 
no output when no photon is detected and a fixed output when one or more photons are 
detected. The action of an on/off detector is described by the following two- value positive 

n 

operator-valued measure (POVM) {no,,,,D, Hi £>} |j| 

^ oo X X A; 

fc=0 ^ 

^l,ri,D = I - ^0,v,D , (lb) 

r] being the quantum efficiency and D the mean number of dark counts, i.e., of clicks with 
vacuum input. In writing Eq. Jl} we have considered a thermal background as the origin 
of dark counts. An analogous expression may be written for a Poissonian background (see 
AppendixEJ. For small values of the mean number D of dark counts (as it generally happens 
at optical frequencies) the two kinds of background are indistinguishable. 

Overall, taking into account the displacement, the measurement on both modes a and b 
is described by the POVM (we are assuming the same quantum efficiency and dark counts 
for both the photodetectors) 

nir(a,/5) = ni'''^)(a)®n^^)(/3) , (2) 

where h,k = 0,1, and Il^i^'^\z) = D{z)Ilh^n,D D^{z), D{z) = exp [za^ — z*a} being the 
displacement operator and z ^ C a complex parameter. 

In order to analyze the nonlocality of the state g, we introduce the following correlation 
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function: 



E,Mc^,f3) = J2 i-)'^' (nir(«>/?)) (3) 

h,k=0 

= 1 + 4 J^,z)(«, P) - 2 [6^,,d(«) + yrjAP)] , 



where 



I,M<y,P) = {^^oo''\<y,P)) (4a) 
g,,o{a) = (u^^'''\a)^l) (4b) 

y,,n{P) = (i^n^^'''\p)) , (4c) 



and where (A) = TT[g A] denotes ensemble average on both the modes. The so-called 
Bell parameter is defined by considering four different values of the complex displacement 
parameters as follows 

Br,,D = ^,,D(a, P) + ^,,D(a', P) + ^,,D(a, P') 

-E„,D{a',P') (5) 
= 2 + 4 {X^,B {a,P)+ Ir,,D («', P) + Iv,D («, P') 

-i,,D{a',p') - g,,D{a) - y.MP)} ■ (6) 

Any local theory implies that satisfies the CHSH version of the Bell me.uality, i.e., 

\^r],D\ < 2 Va, a', P, P' |J|, while quantum mechanical description of the same kind of exper- 
iments does not impose this bound (see Section |IX| for more details on quantum-mechanical 
bounds on in on/off experiments). 

Notice that using Eqs. ^ and ^ we obtain the following scaling properties for the 
functions J^^£i(a,/?), Gn,D{oi) and 3^,,,d(/?) 

^vA'^^P) = (yY^j^ ^v/{^+D)ioi, P) (7a) 

Gv,D (a) = ^v/ii+D) («) (7b) 

3^.,i^(/5) = %/(!+!)) (/3) (7c) 

where = X^ o, Gr] = Gr^fl, and = %,o- Therefore, it will be enough to study the Bell 
parameter for D = 0, namely = i3^,o, and then we can use Eqs. to take into account 



the effects of non negligible dark counts. From now on we will assume D = and suppress 
the explicit dependence on D. Notice that using expression © for the Bell parameter the 
CHSH inequality IBrj^ol < 2 can be rewritten as 

-1 < J^,z)(a,/3) +X^,D(a',/5) + J^,D(a,/5') 

-I^,nia',l3') - g^,nia) - < , (8) 

which represents the CH version of the Bell inequality for our system j^. 

In order to simplify the calculations, throughout this paper we will use the Wigner for- 
malism. The Wigner functions associated with the elements of the POVM (QJ for D = are 
given by (see Appendix EJ 

W[Uo,^]{z) = ^ exp {-A^ \z\^} , (9) 
W[Il,Jiz) = W[I]iz) - W[no,^]iz) , (10) 

with = 2ri/{2 — rj), and iy[I](2;) = vr^^ Then, noticing that for any operator O one has 

W[D{a)OD\a)]{z) = W[0]{z-a), (11) 

it follows that W[D{a) Ho,,, -D''"(a;)](z) is given by 

W[D{a) U,,r,D^{a)]{z) = W[U,,r,]{z - a) , (12) 

and therefore 

W[U^^^'\a,P)]{z,w) = W[Uo,r,]{z - a) 

X W[Ilo,,]{w - (3) (13) 

W[Uo,r,{a) ® I]{z,w) = W[Uo,r,]{z - a) (14) 

W[I^ no,,(/5)] {z, w) = TT-i W[UoJ {w-P). (15) 

Finally, thanks to the trace rule expressed in the phase space of two modes, i.e., 

Tr[Oi O2] = 71^ [ (fz (fw W[Oi] {z, w) [O2] {z, w) , (16) 

one can evaluate the functions X^(a,/?), ^r?(«), and %(/5), and in turn the Bell parameter 
Brj in Eq. as a sum of Gaussian integrals in the complex plane. 
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III. NONLOCALITY OF THE BELL STATES 



We start our analysis by considering balanced superpositions of of zero- and one-photon 
states, i.e., the so-called Bell states, which are described by the density matrices 

P± = |^±)(^'±|, fT± = |<l>±)($±| (17) 

where 

l*±) = ^(|l)|0)±|0)|l)), (18) 

I<^>±) = ^(|0)|0)±|1)|1)). (19) 

In optical implementations Bell states |\l/±) are obtained from single-photon sources using 
linear optical elements, while preparation of requires active devices based on sponta- 
neous parametric down-conversion. 

The Wigner functions of the Bell states are given by 

W[g±]{z,w)= ^ exp{-2|^p -2|w|2} 

X {2\zTwf - 1) , (20) 

and 



W[a±]{z,w)= — exp{-2\z\'^ -2\w\H 

X [1 -2\zTwf + 8\z\'^\w\'^) , (21) 

respectively. 

Let us first consider g±. In this case the functions in Eqs. (jl} are given by 

X.(«, P) = l e-''(l"l'-l^l') [2(1 - r/) + ^ (3\^] (22) 
g,{a) = y,{a) = i e-"!"!' [2 - rj + r/>|2] , (23) 

while the Bell's parameter is obtained using Eq. Maximization of carried out using 
both analytical and numerical methods, indicates that the imaginary parts of the parameters 
a,a',(3,(3' can be neglected for rj = 1, while it influences only slightly the value of \Bn\ for 
T] < 1. More details about the choice of the parametrization are given in Sec. lVIIll Using the 
parameterizations: a = —(3 = J , a' = —(3' = —y/TlJ' for the state |\I'+), and a = (3 = J , 
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a' = = — a/TTj/" for the state \^ -) with JT" G M we get the same Bell's parameter for both 
the states and a maximum violation \B\ = 2.68 (when rj = 1). The Bell's parameter for g± 
is shown in Fig. [2l (a) as a function of J' and r]. 
If we consider a± , we have 



whereas Qn{c() and 3^,,(/3) are given in Eq. (f23|) . As for the states \'^±), the optimal 
parametrization has been obtained by a semi-analytical analysis. We get a = —(3 = J , 
a' = -P' = -y/TlJ for the state |$+), and a = l3 = J,a' = (3' = -y/TlJ for the state 
|$_) (see Sec. IVIIll for more details). Thanks to this choice, is maximized (when r] = 1) 
for both the Bell states a±. The results are shown in Fig. [21(b). 

The overall effect of non-unit quantum efficiency is to reduce the interval of J' values in 
which there is violation. Notice that the states I'I'-i-) are slightly more robust than the \'^±) 
one. In fact, one has \Brj\ < 2. as far as rj falls below 83.6% for |\l/±) and 81.6% for |$±). 
These results are consistent with the study given in Ref. 1^1, where the authors also have 
taken into account mode mismatch and have used a numerical algorithm in order to find the 
best choice of the parameters a, P, a', and /?'. 

A. Unbalanced superpositions 

Our analysis of on-off photodetection is aimed to describe optical implementations of 
nonlocality tests, where most of the experiments have been realized. In this framework the 
Bell states may be obtained from single-photon sources using balanced beam-splitters. 
In order to take into account possible imperfections it is worth to analyze nonlocality prop- 
erties of the class of states that can be obtained from unbalanced beam splitters. Indeed, 
the analysis given above can be extended in order to describe general superpositions of the 
form 



Since the calculations are similar to the ones of the Bell states, here we do not explicitly 
write the analytical results for the states and |$(^). Rather, we plot the corresponding 
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T^(«, f3) =1 e-'^d^l'+l/^l') {2(1 -^)+r]^[l + \a± f]' 




(24) 



= sin¥;|l)|0) + cos(/?|0)|l), 
= sin(^|0)|0) + cos(^|l)|l). 



(25) 
(26) 



Bell parameter Bn in Figs. El and HI In both the plots we used the same parametrization as 
for the Bell states. As one can see in Fig. O in the case of the superposition the best 
result are obtained for the balanced superposition, namely ip = 7r/4. On the other hand, 
the case of |$^) shows a different behavior: here the maximum of the violation for the ideal 
case (i. e. r] = 1) is achieved for a value of ip slightly smaller than 7r/4, and it increases 
as the detection efficiency decreases. Moreover, by the comparison between Fig. [2l (b) and 
Fig. El we can see that, for the particular choice of the parametrization, when rj = 0.8 
the balanced superposition does not violates the CHSH inequality, whereas, adjusting the 
parameter (f, the unbalanced superposition violates it and it does until the efficiency falls 
below the threshold value r] ^ 0.74. 



IV. NONLOCALITY OF SUPERPOSITIONS CONTAINING TWO-PHOTON 
STATES 



A two-photon state which can be easily produced is the one obtained when two single- 
photon states impinge simultaneously in a balanced beam splitter: in this case the output 
state is given by 



|1]) = -^(|2)|0) + |0)|2)) 
and the corresponding Wigner function reads as follows: 

W[n]{z,w)= ^ exp {-2\z\^ - 2\w\^} 



1 /( / I |2 I I |2 12 *2 1 2 

1 — 42; + \w\ — \z — W 



Now, the functions defined in Eqs. (|H) are given by 

[i-v)[i-v + v'{\cy\' + \m] 



X 



(27) 



(28) 



Qr,{a) =%(a) = e' 



-r)|Q!|- 



1-V + 



+ r]^(l-rj)\a\^ + !L\a\' 



(29) 



(30) 
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The parametrization which maximizes the violation of the inequality \Brf\ < 2 is a = /3 = 
and a' = (3'* = ^j2e™l'^J and, as for the Bell states, has been obtained by means of a semi- 
analytical analysis. As it is apparent from Fig. El in this case the violation is quite smaller 
than the previous ones: the maximum violation = —2.07 is achieved when J = 0.45 and 
rj = 1. Also the threshold for violation on the quantum efficiency is higher: for rj < 92% we 
have \B„\ < 2. 



V. NONLOCALITY OF THE TWIN BEAM 



The twin-beam state (TWB) of radiation 

\r) = —^^=^= 2^ tanh" r \n) ® \n) 
V cosh r f—f. 



n=0 

may be produced by spontaneous downconversion in a nonlinear crystal. TWB is described 
by the Wigner function 

4 

WJz, w) = — exp{-2A{\z\'^ + \w\'^) 

+ 2B{zw + z*w*)} , (31) 

with A = A{r) = cosh(2r) and B = B{r) = sinh(2r), r being the so-called squeezing 
parameter of the TWB. Since Wr and the Wigner functions of the POVM ((SJ are Gaussian, 
it is quite simple to evaluate X^(a,/5), ^r?(a), and yriiP) of the correlation function ^ and, 



then, Bri', we have 



4-\4 (r) ~ 
I,{a, (3) = exp {-F, i\a\' + \(3\') 



+ H^{a(3 + a*(3*)} (32) 



with 



F^ = F,(r) = A,-(2A + A^)A^^(r) (33) 
= H,{r) = 2BM,{r) (34) 

^^^'^ = 4iA^-B^H4AA, + Ar ^''^ 
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and 

2A 

g,{a) =y,{a) = ^^^^^^^ 

In order to study Eq. we consider the parametrization a = —[3 = J and a' = 
—(3' = —\f\\J (as in the case of the Bell states, more details are given in Sec. IVIIlll . The 
parametrization was chosen after a semi-analytical analysis and maximizes the violation of 
the Bell's inequality (for r] = \). In Fig.lHlwe plot Br, for r/ = 1: as one can see the inequality 
\Br\ < 2 is violated for a wide range of parameters, and the maximum violation (i3^ = 2.45) 
is achieved when J = 0.16 and r = 0.74. 

The effect of non-unit efficiency in the detection stage is to reduce the the violation; this 
is shown in Fig. where we plot as a function of J' with r = 0.74 for different values 
of the quantum efficiency. Note that though the violation in the ideal case, i.e., r/ = 1, is 
smaller than for the Bell states, the TWBs are more robust when one takes into account 
non-unit quantum efficiency. Comparison between Figs. El and shows that for r] = 0.8 we 
have a region of J' values for which > 2 in the case of the TWB, whereas there is no 
violation for the Bell states. Our parametrization maximize the violation when rj = 1: in 
this way < 2 when r] < 0.77 and r = 0.74. Using different values of a, /?, a', and /3' 
(which, now, depend on rj and the squeezing parameter r), one can extend the violation to 
lower detection efficiency j^. In Sec. IVIIll we will draw some remark about the choice of the 
parametrization. 



VI. EFFECT OF DARK COUNTS 



In the previous Sections we studied the nonlocality of Bell-like states and of the TWB. We 
took into account the quantum efficiency i] and ignored the effects of dark counts: this is a 
quite good approximation, since, at optical frequencies, dark counts may often be neglected. 
However, there are situations in which the effect of dark counts cannot be ignored. In this 
cases, we can add to our analysis the effect of the dark counts using Eqs. ((H). In Fig. IHlwe 
plot Bfj^o for the Bell states and the TWB: as on may expect, the violation is reduced. 

When the number of dark counts is small we can expand the POVMs (jlj and ^ up to 
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first order, arriving at 

n^^) {a) = {l-D- r^Dd,) n^°) (a) (37) 
nS'^)(a,/5) = (1 - 2D - vDd,) n('''°)(a,/?) . (38) 

Now, using Eq. (jT?^ one can express the correlation functions in two equivalent forms as 
follows 

E^,fl(«, (3) = {1-2D- 7]Ddrj) Er,,oia, (5) 

+ 2D [Gr,fi{a) + yr,fl{l3)] (39) 
= {l-D-7]Dd^)E^^o{a,(3) 

-4DJ^,o(a,/3) , (40) 

which, in turn, can be used to express the Bell parameter, as follows 

Br,,D ={l-2D- r]Dd^) 5^,0 (a, P) 

+ iD[gr,,o{a)+y,,om . (41) 

VII. NONLOCALITY OF THE DE-GAUSSIFIED TWIN BEAM 

The de-Gaussification of a TWB can be achieved by subtracting photons from both modes 
0,33. In Ref. Q we referred to this process as to inconclusive photon subtraction (IPS) 
and showed that the resulting state, the IPS state, can be used to enhance the teleportation 
fidelitv of coherent states for a wide range of the experimental parameters. Moreover, in 
Ref. we have shown that, in the absence of any noise during the transmission stage, 
the IPS state has nonlocal correlations larger than those of the TWB irrespective of the IPS 
,ua„tu.efflcie..,eea.oRefs.00). 

First of all we briefly recall the IPS process, whose scheme is sketched in Fig. IHl The 
two modes, a and b, of the TWB are mixed with the vacuum (modes c and d, respectively) 
at two unbalanced beam splitters (BS) with equal transmissivity; the modes c and d are 
then detected by avalanche photodetectors (APDs) with equal efficiency, which can only 
discriminate the presence of radiation from the vacuum: the IPS state is obtained when the 
two detectors jointly click. When the input state, namely the state arriving at the two beam 
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splitters, is the TWB of Eq. (j3T]l . the state produced by the IPS process reads as follows 
(see Ref. [l^ for details) 

where 

P^^(r,T,e) ^±^^^ (43) 
is the probability of a click in both the APDs. In Eqs. (jl^ and we introduced 

C^-C.(r.T..)^-^^^-i^^, (44) 

and defined 

W!:^^^^iz,w) = exp{-Fk\z\^ - Gk\w\^ 

+ Hk{zw + z*w*)} , (45) 

where = {b - fk), Gk = {b - g^), = {2BT + h^), Cu = Ck{e) with Ci = 1, C2 = Cs = 
-2(2 - s)-^, C4 = 4(2 - e)~^; Xk = Xk{r, T, e) and yk = yk{r, T, e) are 

xi = X3 = yi = y2 = a 

X2 = Xi = y3 = y^ = a + 2e{2 - sY^ 

with a = a{r,T) = 2[A{1 -T)+T],b = b{r,T) = 2[AT + (1 - T)]; finally, fk, gu, and hk 
depend on r, T and e and are given by 

fk = Mk [xkB^ 

+ AB\l-A){l-T)+yk{l-Af], (46) 

gk = Uk [xfc(l - Af 

+ AB\l-A){l-T)+ykB^], (47) 
hk=Mk{{xk + yk)B{l- A) 

+ 2B[B^ + Af]{l-T)} , (48) 



The state given in Eq. ()42|1 is no longer a Gaussian state and, in the following, we will use 
the measurement described above in order to test its nonlocality. Nonlocal properties of the 

12 



IPS state ()i2|l have been investigated in Refs. by means of other kinds of nonlocality 

tests. In particular, Ref. addressed the presence of noise during the propagation and 
detection stages, showing that the IPS process onto TWBs is a quite robust method to 
enhance their nonlocal correlations especially in the low energy (i.e., small r) regime. 

In the case of the state the correlation function Q reads (for the sake of simplicity 
we do not write explicitly the dependence on r, T and e) 

4 



k=l 

-2[^f («) + 3^W(/5)]}, (50) 



where 



a\ 



-Ff + («/? + «*/?*) }, (51) 

with ) = Fi''\r,T,e), = Gi^\r,T,e), and H^''^ = H^^''\r,T,e) given by 

) = A, - (Ffc + A,) A<f )(r, T, e) , (52) 

= A, - {Gk + A,) )(r, T, e) , (53) 

= i/fcAlW(r,T,£), (54) 



respectively, and 



"^^P^-G.(F. + A,)-i7|l"l ^ ^''^ 



[Fk (Gfc + A,) - Hi] r, 



xexp^- j^f^-/^^^lJ /?P^ (57) 



(Gfc + A,) - Hi 

In order to investigate the nonlocality of the IPS by means of Eq. (jHI), we choose the same 
parametrization as in Sec. El The results are showed in Figs. [101 and [TT] for 77 = 1 and e = 1: 
we can see that the IPS enhances the violation of the inequality \Br,\ < 2 for small values 
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of r (see also Refs. 0, 12, Moreover, as one may expect, the maximum of violation is 
achieved as T ^ 1, whereas decreasing the effective transmission of the IPS process, one 
has that the inequality becomes satisfied for all the values of r, as we can see in Fig. [TT] for 
T = 0.8. 

In Fig. ini we plot for the IPS with T = 0.9999, e = 1 and different r]. As for the 
TWB, we can have violation of the Bell's inequality also for detection efficiencies near to 
80%. As for the Bell states and the TWB, a rj- and r-dependent choice of the parameters 
in Eq. © can improve this result. The effect on a non-unit e is studied in Fig. where 
we plot Bri as a function of T and e and fixed values of the other involved parameters. We 
can see that the main effect on the Bell parameter is due to the transmissivity T. 

The presence of dark counts at the detection stage can be taken into account using 
Eqs. 0: since the results are similar to those of the Bells states and the TWB presented in 
Sec. IVIl we do not report them explicitly. 



VIII. CHOICE OF THE PARAMETRIZATION 



In this Section we draw some remark about the choice of the parametrization used in 
the investigation of the Bell parameter Br,. Numerical analysis has shown that, in the case 
of state |\E'+) and in the presence of non-unit quantum efficiency, the maximal violation of 
the Bell's inequality for the displaced on/off test is achieved by choosing a, a', (3, and /3' 
as complex parameters 1^. On the other hand, here we addressed only real parametrization 
for the Bell's parameter B given in Eq. (jH)), and, in particular, we take a = —[3 = J7 G M, 
and a' = -/?' = J^' g M for the states |^+), 1$+), the TWB and the IPS state, while we 
put a = /? = G M and a' = = J^' G M for the states |^'_) and |<l>_). In Fig. d we plot 
\Bn\ as a function of J and J' in the ideal case (i.e., r] = 1) for (a) the states \^±) and (b) 
the TWB with r = 0.74, which maximizes the violation. The results for the other states are 
similar. In both the plots, the darker is the region, the bigger is the violation (the white 
region refers to < 2). As one can see, there is a symmetry with respect to the origin, 
which implies that the best parametrization has the form J' = —kJ^ with k G M, k > 0. 
Furthermore, for all the considered states, the numerical analysis shows that a good choice 
for K is K = vTT, which is an approximation of the actual value. 

In Fig. [ini the effect of ri is taken into account: Since the results for the Bell's states, 
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the TWB and the IPS state are similar, we only address the TWB case: there is still the 
symmetry with respect to the origin, but a thorough numerical investigation shows that the 
maximum of Br,, and, then, k depend on both r] and r. 

Notice that we have considered real values for the parameters. It can be shown numeri- 
cally that for decreasing rj a complex parametrization leads to a slight improvement. 



IX. BELL'S INEQUALITY, POVMS AND MAXIMUM VIOLATION 

In this Section we address the maximal violation of the Bell inequality that is achievable 
by using non-projective measurements. 

Let us consider two systems, A and B, and the generic POVM {no(C), ni(C)}, depending 
on the complex parameter (, such that ^i{C) = I — no(C)- We define the observables 

^k(C) = ni(C) - no(C) = 1-2 no(C) (58) 

acting on system k = A, B, respectively (we are using the same POVM for both the systems), 
i^'urthermore, we assume that ^a(C) and ^r(C) have spectra included in the interval [—1, +1] 
1^ . Now we introduce the Bell operator 4| 

^ =^a(«) ® MP) + ^a(«') ® MP) 

+ ^A(a) ® W) - ^A(a') ® W) , (59) 

which has the property [3| 

v2 

+ ^l(a')®I + I® ^b(/3')] ■ (60) 

If ni(C) and no(C) are projectors on orthogonal subspaces, namely n^(C) = ni(C), ^o{() = 
no(C) and no(C) ni(C) = 0, then 

^a(C) = ^b(0 = I, (61) 

and Eq. (|6n|) leads to 

B<2V2, (62) 

where B = Ti[q^], g being the state of the system, is the Bell parameter. The bound 2a/2 
is usually known as Cirel'son bound and is the maximum violation achievable in the case of 
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a bipartite quantum system Eq. (|62|) may be also derived in a different since 
tfie squared Bell operator reads 



S2 



41 + [^A(a), ^a(«')] ® [^b(/3), MP')] , (63) 



then using the relation ||[v4, < 2||y4|| \\B\\, where ||y4|| = Sup|^^||y4|V^)||, we have B'^ < 8, 
from which Eq. |62|l follows. 

On the other hand, when {no(C), ni(^)} is not a projective measurement a different 
inequality should be derived. First of all we note that the observables ^k(C) corresponding 
to the POVM given in Eq. ^ satisfy the hypothesis of the Cirel'son theorem. In fact, in 
this case 

oo 

^k(C) = 5^ [1 - 2(1 - r^Y] D{C)\n){n\D\0 (64) 

n=0 

and its spectrum {A^}, = 1 — 2(1 — 77)", lies in the interval [— 1, +1] for < < 1 (when 
T] = 1 the spectrum reduces to the two points {—1, +1})- Now, one has 

ffl{0 = 1-4 no(C) ni(c) = I - 4 [no(C) - n^(c)] 

= I-4^a(C), (65) 
and, analogously, ^b(0 = I — 4<fB(0i where we defined the operator 

4(C) = no(C)-n^(c), (66) 

k = A, B. In this way, from Eq. (|6n|) follows 

^ < 272(1 - [^A(a) + + '^a(«') + W)]] , (67) 

with A(C) = A(C) ® I and (^6(0 = I ® ^^6(0- Finally we get 

i3<i3('^^")(a,/3,a',/3'), (68) 

where we defined 

= 2v^{1-[^a(«)+^b(/5) 

+ ^a(«)+^b(/30]} (70) 

and £a(C) = Tr[<fA(C) ® I] and ^b(0 = Tr[I ® ^b(0]- Now, since 4 > 0, k = A, B, one has 
that the bound of the Bell parameter is smaller than the limit 2v^, obtained in the case of 
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projective measurements. Notice that the new bound depends on the parameters a, (3, a', 
and 13' of the measurement and on the state under investigation itself. 

In the following we address the problem of evaluating the maximum violation for the Bell 
states, the TWB and the IPS state when a non-unit efficiency affects the displaced on/off 
photodetection. First of all we note that 

QO 

nM = E(i-^)''l^)<^l=no,.(2-.), (71) 

A:=0 

SO that 

£K{0 = Q,{C)-Gv{2~v)iO (72) 

^b(0 = %(C)- 3^.(2-.) (0- (73) 

In this way it is straightforward to evaluate B^'^'^^ for the Bell states, the TWB and the IPS 
state. The results are shown in the Figs.HBHSU in Figs.HHHni and ED we plot bI^^^^ using 
the parametrization a = (3 = J and a' = (3' = vTTj", which maximizes the Bell parameter 
Br,. It is worth noticing that for all the considered states and for fixed J the limit b'^'^^^ 
is never reached; on the other hand, even if the actual maximum violation, i.e., i3^,max, is 
quite lower than the Cirel'son bound 2\/2, it is relatively near to the new bound given by 
Eq. |68|l . 

Notice that a similar analysis may be performed through the squared Bell operator, which 
for ^k(C) 7^ I is given by 

-[^A(a),^A(«')]+® (^b(/?)-^b(/5')) 

+ [^A(a), ^a(«')] ® Wb{P) + ^b(/3')] , (74) 

[A, = AB + BA being the anti-commutator. As for Eq. (|68| . the maximum value of 
Eq. (f7i|l depends on the state under investigation and on the POVM itself. 

X. CONCLUDING REMARKS 

We have analyzed in details the nonlocality of several two-mode (entangled) states of light 
by using a test based on displaced on/off photodetection. Nonlocality has been quantified 
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through violation of CHSH inequality for the Bell's parameter. Effects due to non-unit 
quantum efficiency and nonzero dark counts have been taken into account. We found that 
unbalanced superpositions show larger nonlocality than balanced one when noise affects 
the photodetection process, and that twin-beam nonlocality is more robust than that of 
superpositions of few photon-number states. De-Gaussification by means of (inconclusive) 
photon subtraction is shown to enhance nonlocality of twin beams in the low energy regime. 
We have also shown that, since our measurement is described by a POVM rather than a 
set of projectors, the maximum violation the CHSH inequality cannot saturate the Cirel'son 
bound. A novel state-dependent bound has been derived. 

Acknowledgments 

Fruitful discussions with A. Ferraro are kindly acknowledged. This work has been par- 
tially supported by MIUR (FIRB RBAU014CLC-002). 

APPENDIX A: NOISY ON/OFF PHOTODETECTION 

The action of an on/off detector in the ideal case is described by the two- value POVM 
{ITo = |0)(0|, III = I — Flo}, which represents a partition of the Hilbert space of the signal. In 
the realistic case the performances of on/off photodetectors are degraded by two effects. On 
one hand, one has non-unit quantum efficiency, i.e., the loss of a portions of the incoming 
photons, and, on the other hand, there is also the presence of dark-count, i.e., by "clicks" 
that do not correspond to any incoming photon. In order to take into account both these 
effects we use a simple scheme described in the following. 

A real photodetector is modeled as an ideal photodetector (unit quantum efficiency, no 
dark-count) preceded by a beam splitter (of transmissivity equal to the quantum efficiency 
77) whose second port is in an auxiliary excited state which can be a thermal state, or 
a phase-averaged coherent state, depending on the kind of background noise (thermal or 
Poissonian). If the second port of the beam splitter is the vacuum u — |0)(0| we have no 
dark-count; for the second port of the BS excited in a generic mixture u — the 
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\n){n\ . 



(Al) 



POVM for the on/off photodetection is given by (Hi = I — Ho) 

OC OO y 

n. = E(i-'')"E''..''' ( " Vl 

The density matrices of a thermal state and a phase-averaged coherent state (with M mean 
photons) are given by 

' M ' 



M + 



+ 1^ 



M+ 1 



\s){s\ 



Up — e 



(A2) 
(A3) 



In order to reproduce a background noise with mean photon number D we consider the state 
u with average photon number M — D/{1 — rf). 
In this case we have 



\n) {n\ 



^l,r,,D - e ^E 



1 — 1] 



\n){n\ , 



(A4) 
(A5) 



where t and p denotes thermal and Poissonian respectively, and Ln{x) is the Laguerre poly- 
nomial of order n. The corresponding Wigner functions are given by 



n 2(1 + D)-r] 



TT 2 — T] 



exp 



exp 



_2iT_ 
2-77 
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2(1 + D)-ri 
{D + \a\')] /o 



\a\ 



4:\a\y/rjD 



(A6) 
(A7) 



respectively, where Io{x) is the 0-th modified Bessel function of the first kind. For small D 
the POVMs coincide up to first order, as well as the corresponding Wigner functions. 
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FIG. 1: Scheme of the nonlocality test based on displaced on/off photodetection: the two modes a 
and 6 of a bipartite state g are locally displaced by an amount a and /? respectively, and then revealed 
through on/off photodetection. The corresponding correlation function violates Bell's inequalities 
for dichotomous measurements for a suitable choice of the parameters a and /?, depending on the 
kind of state under investigation. The violation holds also for non-unit quantum efficiency and 
non-zero dark counts. 
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FIG. 2: Plot of —Brj for the states \^±) (a) and of Br^ for |<I>±) (b) as functions of JT" and rj. The 
maximum violations for r/ = 1 are: (a) —Brj = 2.69, and (b) Brj = 2.68, which are both obtained 
when J' = 0.17. For the particular choice of the parametrizations Br^ is the same for |^'±) and for 
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(b) ip/TT 

FIG. 3: Plot of —Brj: (a) for |^^) as a function of J and ip in the case of ideal (i.e., r] = 1) on/off 
photodetection; (b) for J = 0.17 and different values of rj: from bottom to top rj = 1.0, 0.9, 0.85, 
and 0.8. The vertical line is (f = it /A. 
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(b) ^p/tt 



FIG. 4: Plot of Br^: (a) for \^ip) as a function of J and (p in the case of ideal (i.e., rj = 1) on/off 
photodetection; (b) for J = 0.17 and different values of r/: from top to bottom rj = 1.0, 0.9, 0.85, 
0.8, and 0.75. The vertical line \s ip = 7r/4. 




FIG. 5: Plot of —Br^ for the superposition of two photons as a function of J and r]. The maximum 
violation is —Br^ = 2.07, which is obtained when J = 0.45 and rj = 1. 
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FIG. 6: Plot of Br, for a TWB as a function of JT" and the TWB squeezing parameter r in the case 
of ideal (i.e., rj = 1) on/off photodetection. The maximum violation is Br, = 2.45, which is obtained 
when J = 0.16 and r = 0.74. 
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FIG. 7: Plot of Br, for a TWB as a function of ^7 with r = 0.74 for different values of ij: from top 
to bottom r) = 1.0, 0.9, 0.85, and 0.80. 
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FIG. 8: Plot ofBrj^oia, /?, /?') for the Bell states |$±) (upper plot) and the TWB (lower plot). We 
set J' = 0.17 for the Bell states, J' = 0.16 and r = 0.74 for the TWB, and used the parametrizations 
introduced in Sees. IIIII andlv ] respectively. 




FIG. 9: Scheme of the IPS process. 
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FIG. 10: Plot of Brj for the IPS state with T = 0.9999 and e = 1 as a function of J and the 
TWB squeezing parameter r in the case of ideal (i.e., ?] = 1) on/off photodetection. The maximum 
violation is Brj = 2.53, which is obtained when J7 = 0.16 and r = 0.39. 
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FIG. 11: Plot of for the IPS state as a function of J7 with r = 0.39 for different values of T and 
£ = 1 in the ideal case (i.e., rj = 1): from top to bottom T = 0.9999, 0.95, 0.90, and 0.80. 
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FIG. 12: Plot of Br, for the IPS state as a function of J with r = 0.39, T = 0.9999, £ = 1, and for 
different values of r]: from top to bottom rj = 1.0, 0.9, 0.85, and 0.8. 
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FIG. 13: Plot of for the IPS state as a function T and e with J = 0.16, r = 0.39, and, from top 
to bottom, r] = 0.99, and 0.90. The main effect on is due to the transmissivity T. 
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FIG. 14: Plots of \Br,\ for (a) the states |*±) and (b) the TWB with r = 0.74 and r) = I. The 
darker is the region, the bigger is the violation of the Bell's inequality. In the white region < 2. 
J7 and J'' refer to the particular parametrization of B^j, see the text for details. 
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FIG. 15: Plots of in the case of TWB for different 77. The TWB parameter r is chosen in 
order to maximize the violation of the Bell's inequality as 77 varies. We put (a) = 1, r = 0.74 
(B^,ma. = 2.45), (b) T] = 0.9, r = 0.67 {Br,,ma. = 2.24), (c) rj = 0.85, r = 0.60 (B^,^ax = 2.15), and 
(d) r] = 0.8, r = 0.49 (6,j,max = 2.07). J' and J'' refer to the particular parametrization of Br^, see 
the text for details. Since there is symmetry with respect to the origin, we show only the region 
J' > 0. 
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FIG. 16: Plot of B)j"^>{a,P,a',P') for the Bell states. We set a = p = J and = = J' and: 
(a) T] = 0.9, (b) r] = 0.8. 
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FIG. 17: Plots oiBi'^^\a,(3,a',(3') in the case of the Bell states for different values of rj; from top 
to bottom (solid lines): rj = 0.95,0.9,0.85, and 0.8. The dashed line corresponds to the value 2\/2 
obtained when r] = 1. We set a = /3 = JT" and a' = (5' = y/VVJ, which maximize B. 
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FIG. 18: Plot ofBlj"^' {a, (3, a', /?') for the TWB with r = 0.74. We set a = ^ = J and a' = (3' = J' 
and: (a) ry = 0.9, (b) r] = 0.8. 
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FIG. 19: Plots of bJ™^^ (a, ^, a', in the case of the TWB with r = 0.74 for different values of r?; 
from top to bottom (solid lines): rj = 0.95,0.9,0.85, and 0.8. The dashed line corresponds to the 
value 2\/2 obtained when r/ = 1. We set a = (5 = J and a' = f3' = \p\\J^ which maximize B. 
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FIG. 20: Plot ofBi^^'{a,P,a',P') for the IPS state with r = 0.39, T = 0.9999, and £ = 1. We set 
a = P = J and a' = P' = J' and: (a) r] = 0.9, (b) ry = 0.8. 
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FIG. 21: Plots of fij"'"''' (a, /?, a', in the case of the IPS state with r = 0.39, T = 0.9999, 
and e = 1 for different values of r/; from top to bottom (solid lines): ry = 0.95,0.9,0.85, and 0.8. 
The dashed line corresponds to the value 2\/2 obtained when rj = 1. We set a = /3 = J7 and 
a! = (3' = \/lTi7, which maximize B. 
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